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Abstract In this work we apply the formalism developed in [M. Lepers et al, Phys. Rev. A 77, 043628 
(2008)] to different initial conditions corresponding to systems usually met in real-life experiments, and 
calculate the observable quantities that can be used to characterize the dynamics of the system. The 
position space point of view allows highly intuitive pictures of the physics at play. 
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1 Introduction 

Over the last 30 years, the kicked rotor has played a cen- 
tral role in the field of quantum chaos, which is defined as 

. the study of quantum systems whose classical counterpart 

: is chaotic [I][2J[3]. The first atom-optics experimental re- 
alization of its quantum version in 1995 [3] was obtained 
by submitting laser-cooled atoms to a periodic series of 
"kicks" of a laser standing wave which is put on periodi- 
cally for a very short time compared to the atom dynam- 
ics. It has been shown that, under quite general conditions, 
this system displays a characteristic behavior named "dy- 
namical localization", that is, an asymptotic saturation of 

\ its average kinetic energy, that can be attributed to the 
existence of destructive interferences in the system. An im- 
pressive number of experimental results followed this pi- 

; oneering work [5l fBlf7lfg|J^lfTUlfTTlfT^lfT5] . spanning domains 
like quantum transport |14| . measurements of the grav- 
itational constant [15] . condensed matter and quantum 
phase transitions with, recently, the first experimental ob- 
servation of the metal-insulator Anderson transition with 
matter waves [T31ITtjl[l"T] . 

Dynamical localization is not the only interesting phe- 
nomenon displayed by the kicked rotor. Another strik- 
ing feature of the quantum kicked rotor (QKR) is the 
existence of quantum resonances (QRsjJ which arise if 
the kicking period matches the evolution of the quantum 
phase of the atoms. A QR is characterized by a ballistic 



a Present address: Laboratoire Aime Cotton, Universite 
Paris-Sud, Bat. 505, Campus d'Orsay, F-91405 Orsay Cedex, 
France 

1 Not to be confused with classical resonances or acceler- 
ator modes. These are observed for particular values of the 
stochasticity parameter (see below) K « lira and of the initial 
conditions. 



behavior in which the velocity increases linearly and the 
average kinetic energy quadratically with time, in sharp 
contrast with dynamical localization. Since the pioneering 
theoretical works [T81I19) . the QRs have been the subject 
of an impressive amount of work, both theoretical [20 , 21 , 
[2~2l^3,24,25.26.27'P51[2'9l and experimental [4~PDll3T] [31 
33,34,35,36,37,38], including the effect of gravity (the so- 
called quantum accelerator modes) [7,39,40,41,42,43,44], 
directed motion (or quantum ratchets) [45U46U471I481I14] . 
or the effects of interatomic interactions in a Bose- Einstein 
condensate (BEC) gHEE]. 

The main goal of activity around QRs is to calculate 
and measure the average momentum and kinetic energy, 
and possibly control the dynamics, e.g. by constructing 
ad hoc position or momentum distributions. The present 
article is placed in this mainstream. We apply analytical 
methods previously developed in |27j . validated by com- 
parison with numerical simulations, to the mean momen- 
tum and kinetic energy of the QKR in situations often 
encountered in the laboratory. Some of these results have 
already been presented in other works, but our approach 
and the intuitive picture of the dynamics sheds a new light 
on these complex behaviors. 

Systems that are periodic in space, as the kicked rotor 
is, conserve quasimomentum. This means that its eigen- 
states are combs displaying the same periodicity as the po- 
tential, both in the position and in the momentum space, 
the Fourier transform of a comb being a comb. A comb- 
shaped spatial wave function, in the absence of a potential, 
has the striking property of resuming its initial (strongly 
localized) shape after some characteristic time, an effect 
known in optics as the Talbot effect [51] . The QRs of the 
kicked rotor are a quantum version of the Talbot effect. In 
[2~T] we performed an analysis of such effect in the position 
space, which allowed us to simply relate the QR dynamics 
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to local properties of the potential at the positions where 
the wave function reconstructs. The QRs' ballistic behav- 
ior then arises as a consequence of the fact that the wave 
function focuses, before each kick, at the same position 
with respect to the potential. The system thus receives 
the same amount of momentum from each kick, and the 
effect of successive kicks adds constructively. The result- 
ing dynamics, although based on intereferences, can be 
interpreted in terms of purely mechanical arguments. 

However, the initial conditions of real experiments are 
seldom comb-shaped functions. The initial state of a Bose- 
Einstein condensate (BEC), for example, can be modeled 
as a sharp superposition of plane waves around zero mo- 
mentum. The hotter atom cloud from a magneto-optical 
trap is represented as an incoherent mixture of plane waves. 
The purpose of the present article is to show that the pic- 
ture presented in [27] is also valid in those situations. The 
main point is that an arbitrary initial wave function can be 
viewed as a superposition of comb-shaped states, which, in 
the specific case of QRs, evolve independently from each 
other. We consider in the present work only the case of 
the so-called simple quantum resonances, in which the re- 
constructed wave packet is identical to the initial one. The 
more complex case of the so-called "higher-order" quantum 
resonances, in which the initial wave packet reconstructs 
into many replicas, will be discussed in a future work. Thus 
in what follows, unless otherwise stated, it is understood 
that the term quantum resonance means simple quantum 
resonance. 

The article is organized as follows. In section [3J we 
recall the essential features of our approach to QRs in 
position space. We then generalize these results to initial 
conditions of experimental interest: After considering the 
case of a plnae wave (sec. [3]), which is the basic component 
to describe experimental situations, we address the case 
of a coherent superposition of two plane waves (sec. 2]), 
which is used to interpret the recent experiments on quan- 
tum transport and ratchets, the case of a narrow momen- 
tum distribution, which is characteristic of a BEC (sec. 
[S]); finally, for the sake of completeness, we briefly discuss 
the "trivial" case of a broad initial momentum distribu- 
tion, characteristic of a thermal cloud (sec. [5]). In all these 
cases, our picture of the evolution seen in position space 
allows a simple interpretation of the dynamics. Section [7] 
draws the conclusions of this work. 



2 Quantum resonances in position space 



The atomic kicked rotor is obtained by placing laser-cooled 
atoms of mass M and momentum p in a standing wave 
formed by two counterpropagating (along the x-axis) laser 
beams of wave number — 2tt/\l. The atoms are thus 
submitted to a mechanical potential called "optical poten- 
tial" V(x) = Vq cos (2kzx) formed by the standing wave. 
If the standing wave is modulated in the form of periodic 
pulses (period T) of duration r very short at the timescale 



of the atom dynamic^], one obtains the well-known kicked 
rotor Hamiltonian 

H = — + KcosX^25(t-n), (I) 

~ n 

where we used reduced variables [4] in which time is mea- 
sured in units of the kicking period T, X = 2kL%, P = 
kp/2hk L B K = W T/h, where k = Ahk 2 L T/M plays the 
role of a reduced Planck constant 0. The parameter k. 
which plays a critical role in what follows, can be experi- 
mentally tuned by changing the kicking period T. 
The quasimomentum /3 is defined by the relation 

P=(n + (3)k, (2) 

where n is an integer and f3 € [— ^; |) is in the first Bril- 
louin zone. The spatial periodicity of the potential implies 
that (3 is a constant of motior0. In position space, wave 
functions ij}f)(X) corresponding to a well-defined quasimo- 
mentum have the well-known Bloch-wave structure 

MX) = e^ x u (X), (3) 

where up(X) has the same spatial period 27r as the poten- 
tial 0. 

The dynamical evolution of the atoms is governed by 
the one- (temporal) period evolution operator of the quan- 
tum kicked rotor 

U = cxp (-i^ cos X) CX P ' ( 4 ) 

where the first term on the right-hand side corresponds 
to the effect of an instantaneous kick and the second to 
a free propagation between two kicks. The factorization 

2 That is pr/M <C \l for all relevant values of the atomic 
momentum p. 

3 As the atoms move on a straight line, X G] — oo; +oo[ and 
P e] - oo; +oo [. 

4 In the sense that the corresponding Schrodinger equation 
is ikdip/dt = Hip. 

5 The conservation of quasimomentum is a general property 
of periodic potentials, but the case of the atomic kicked rotor 
allows a completely different interpretation. Interacting with a 
far-detuned standing wave, an atom can only absorb and re- 
emit photons by stimulated emission. Four kinds of process are 
possible: Two of them correspond to the absorption and emis- 
sion of a photon in the same propagative wave forming the 
standing wave; in this case the overall momentum exchanged 
between the atom and field is zero. Two other process corre- 
spond to the absorption of a photon in one of the propagative 
waves followed by the emission of a photon in the other prop- 
agative wave, in which case the momentum exchanged is two 
times the photon momentum, that is 2hkL- Starting from a 
well-defined momentum po only momentum states of the form 
2nhkL (or nk in our units) are populated. The conservation of 
quasimomentum can thus be deduced from this "microscopic" 
argument, not directly related to the periodicity of the poten- 
tial. 

6 The atom wave function is thus ift(X) = J_y 2 df3ipp(X) 
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is allowed because the kicks are short compared to the 
atom dynamics, so that one can neglect the contribution 
of the free-evolution term during this short time. (General) 
quantum resonances appear if k is a rational multiple of 
47r, i.e. k = Airr/s. In physical units, this corresponds to 
a kicking period T = rTx/s, where Tt — irAI/hk^ is a 
characteristic time, called the Talbot time. Given a state 
tpp(X,t) of well-defined quasimomentum /3, one can use 
Eq. ((3]) to calculate the wave function ipp(X,t + 1) one 
period later. The condition for simple quantum resonances 
is k = 2ir£, with t integer, which leads to a particularly 
simple recursion relation [see [57], Eq. (10)] for the wave 
function: 



1>p(X,t) = e -' KcosX / k exp {ikf3{fi + l)/2) 
x x/jp (X~k((3 + l/2),t-l). 



(5) 



The term exp (ikfi{f3 + l)/2) in Eq. © is a global phase 
that does not affect the dynamics. The free evolution be- 
tween kicks thus preserves the shape of ipp but translates 
it, over a whole kick period, of 



(6) 



In [57] we have considered the evolution of an atom 
initially prepared in a Bloch wave perfectly localized in 
a spatial period [— tt, tt) Q and centered at a position £ G 
[— 7r, 7r), that is 



Xp(X-Z) = e iflx £ 6(X-S-2nj) 



(7) 



In what follows, X/3i which will play a central role in our 
interpretation of QRs, will be referred to as a comb-shaped 
Bloch wave (CSBW). In [57| we have shown that the mean 
position and momentum evolution of a CSBW obeys a 
simple map 



^/j,e(*) = - p w(*-i)- 



Ksm^ + vpt) 



+ vps), 



(8) 



(9) 



and the evolution of the mean kinetic energy turns out to 
be given by 



1 



Ep^t) - Ep i6 (0) = - (Pl ( (t) - P^(0)) 



(10) 



Unlike the e-classical map, used to describe quasi-resonant 
dynamics [3"9"1H2"] . Eqs. (|5|)- (flU|) do not describe a classical 
motion, as the change in position Xp^(t) — Xp^(t — 1) 
does not depend on momentum at kick t — 1. 



7 Let us stress that here "localized" means localized in the 
"spatial" first Brillouin zone [— tt,tv), the full wave function is 
a comb of delta functions. In the momentum space the wave 
function is also a Dirac comb; the Heisenberg's uncertainty 
relations are thus fully satisfied. 



The CSBW in Eq. (J7J) is characterized by its quasimo- 
mentum P and by its "center" £ which is determined by 
the weights of each momentum component nk. In the case 
of a quantum resonance - and only in this case - the evo- 
lution between two subsequent kicks just translates the 
wave packet of Vp\ therefore, for fixed j3, contributions 
due to different initial values of £ are not mixed, and each 
CSBW Xp(Q evolves independently of the others. An ar- 
bitrary initial state can be expanded in CSBWs xp(X — £) 
[see Eq. ©] 



Mx,t = o) 



with 



(12) 



ipp(n) being the amplitude of the component n in momen- 
tum space. 

In the following sections, we apply the above formalism 
to cases of experimental interest where the initial state 
is not a CSBW, but either a coherent or an incoherent 
mixture of CSBWs. As the CSBWs evolve independently, 
we have just to find the weights irp £ of each CSBW in the 
initial state, which then allows us to directly calculate the 



averages 



1/2 



P k (t) - P k {t = 0) = / d/3 

J -1/2 J -TT 

x(P^(t)-P* 4 (t = 0)) (13) 

with k — 1,2, and P% + given by Eqs. © and ([IU]). In the 
case where the initial state is a coherent superposition of 
plane waves, np,£ is given by 



= l^(£,* = o)| 2 



while in the case of an incoherent mixture, irp^ is 



(14) 



(15) 



3 Plane wave 

Let us first apply the previous results to the simple case 
of a plane wave. A plane wave of momentum po has a 
well-defined quasimomentum /3q which satisfies po = (tiq + 
f3o)k, but is completely delocalized in the real space, which 
means that all values of £ are equiprobable on [— 7r, tt), 
hence 7173,5 = S((3 - ^)/2tt. From Eqs. © and Q3| one 
finds. 
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Figure 1. (color on line) Dynamics of a zero-momentum plane 
wave. The initial distribution, uniform in the position space 
(restricted here to two potential wells), is divided into "slices" 
represented in different colors. Panel (a) corresponds to the 
ballistic evolution of the quantum resonance k = 4tt. After a 
free evolution between kicks, each color has come back to the 
initial position (slices of same color are at the same position in 
the bottom "degrade" as in the top one). The subsequent kick 
is thus applied to the same spatial distribution, and each slice 
receives the same momentum from any kick. The momentum 
transfers thus add constructively from kick to kick. Panel (b) 
corresponds to the "anti-resonance" condition £ = 2n. In such 
case, each colored slice is shifted by half the lattice step by the 
free evolution between kicks, and thus receives from a given 
kick a momentum transfer of the same magnitude, but of op- 
posite direction with respect to the previous kick, which leads 
to a periodic behavior. 



The average kinetic energy is obtained in the same way: 



function receives the same momentum transfer K sin £ [for 
example, the yellow part receives the maximum momen- 
tum transfer +K, see Eq. ©]. As the wave function fills 
uniformly in the potential wells when the kicks are applied, 
the average momentum is always zero. But, the average 
kinetic energy grows quadratically, since some "colors" are 
accelerated in the positive direction, while others are ac- 
celerated in the negative direction. Although the position- 
space distribution is uniform, the local effect of the kicks 
remains visible on the average kinetic energy. By contrast, 
figure[T](b) illustrates the so-called anti-resonant case, for 
which the position-space kicks, a given "color" receives a 
momentum equal to —Ksin^, whereas at even kicks, it 
receives -t-ifsm^, hence a periodic behavior. 



4 Superposition of two plane waves: 
inducing directed motion 

It has recently been show that directed motion of the wave 
packet (which is sometimes called a "quantum ratchet") 
can be obtained by preparing a BEC in a coherent super- 
position of two momentum states via a Bragg pulse |45[ 
|4"7I|48] . This yields a superposition of two plane waves: say, 
one centered at P — and one centered at P — — k, each 
plane wave being equally populated. We present here an 
alternative, and hopefully more intuitive, interpretation of 
this result using the position-space point of view. 

A coherent superposition of two plane waves (|0) 



ic 



-k) )/ \[2 is characterized by a mean momentum equal 
to -k/2 [IB] . The corresponding initial CSBW distribu- 
tion np^ is obtained from 



V2 



ie 



(18) 



E(t) = E(0) 



47T 



k 2 (n + /3 ) 



\s=l 

K 2 sin 2 (^ t/2) 
4 sin 2 (^ /2) 



(17) 



a result previously obtained in [23J. If vp — 0[2n], where 
x[2ir) means "x modulus 27r", the dynamics is ballistic: 
E(t) = E(0) + K 2 t 2 /4:, whereas if v Pa = tt[2tt] (e.g. k = 
2tt,Pq = 0) the dynamics is periodic (or "anti- resonant") , 
i.e. E(t) = E(0) + K 2 /4 for t odd and E(t) = E{0) for t 
even. 

Figure [1] gives a simple interpretation of this behav- 
ior for a zero-momentum plane wave (no — /3q = 0). On 
panels (a) and (b), the initial position-space distribution 
is uniform, but we have attributed a different color to 
each slice of it, so that one can follow its evolution. Eq. 
([5]) shows that the free evolution between kicks results in a 
simple global shift of the wave function between t = and 
t = 1 kick. In the ballistic case (a), this shift is such that 
all "colors" (or slices) come back to the same place at t = 1 
(and also at t = 2, 3, ...). Therefore, each part of the wave 



and Eq. (jl2j) : 



^, ? H^(^ = o)| 2 = ^(i 



sin(£-0)), (19) 



which is localized in quasimomentum and inhomogeneous 
in the position representation. Performing the integral in 
Eq. COl we retrieve Eq. (5) of ref. gBJ 



. k Kt 
P(t) = ---_ cos < 



(20) 



which shows a directed motion, which is not due to the 
non-zero initial average momentum, but to the inhomo- 
geneity of the position-space distribution. This point can 
be stressed by considering an initial superposition of two 
arbitrary plane waves (| nk) — ie % ^ \(n — 1) £)) j\p2. By ap- 
plying Eqs. (fT4"]) and (fl"2"|) . one obtains the same distribu- 
tion irp^ as in Eq. (fTT)|) . For all value of n, the directed mo- 
tion will be characterized by the same current, whatever 
the initial average momentum and kinetic energy. This re- 
markable property emphasizes the central role of irp^ on 
the dynamical evolution of the QKR. 
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t (kicks) 

Figure 2. Time evolution of the average kinetic energy ob- 
tained by integration of Eq. Q. The initial momentum distri- 
butions are coherent Gaussian superpositions of plane waves, 
whose central quasimomenta j3 — are ballistic (k = 4ir). Their 
root mean squares, a — 0.0115 (solid lines) and a — 0.00866 
(dashed lines), are such that A = 0.04 and 0.03 respectively 
(see text). Other parameter: K = 10. 

In Eq. ([2"0")l , the direction of the motion is controlled by 
the quantum phase. For = 0, the maximum of tt^ is lo- 
cated at £ = —tt/2. where the local momentum transferred 
from the pulse reaches its minimal value — K (cf. Fig. [I]), 
which creates a dominant current towards —X, thus corre- 
sponding to a negative average momentum. The opposite 
situation = tt corresponds to the case where the maxi- 
mum of 7Ty3^ coincides with the maximum positive trans- 
ferred momentum (+K): the current is maximum towards 
the +X direction. The prefactor 1/2 in the current comes 
from the averaging over the width of the spatial distri- 
bution, which means that some slices of the distribution 
receive an amount of momentum which is smaller than 
K in absolute value. The amplitude of the current can 
thus be increased by building a position-space distribu- 
tion sharply localized on favorable values of the spatial 
potential. 

As pointed out in Ref. [45], such kind of quantum 
transport cannot be observed if the initial state is a plane 
wave. Our picture of QRs makes this result easily under- 
standable, as a plane wave equally samples all parts of 
the kicking potential. This reasoning is straightforwardly 
generalized to the arbitrary kicking potentials obtained by 
adding higher harmonics of the spatial frequency 2k l [52] . 

In the anti- resonant case k = 2tt£ {I odd), no transport 
can be observed, even with an inhomogeneous distribu- 
tion, as given slice at position £ will always receive from 
kick t a momentum exactly opposite to the one received 
at kick t — 1 . 



5 Narrow initial momentum distribution 

The momentum distribution of a BEC, although rather 
sharp, has actually a finite width, which induces a damp- 
ing with time of the effects observed in the ideal, zero- 
width case [H]. This is illustrated on Figs. [2] and with 
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3 0.1 
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t (kicks) 
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Figure 3. Same as Fig. [5J except that the central compo- 
nent of the initial momentum distributions is anti-resonant 
(k — 2n,/3 — 0). Their root mean squares are: (a) a = 0.0115 
and (b) a = 0.00577. They are such that A = 0.04 and 0.02 
respectively (see text). 

numerical simulations made for a coherent initial momen- 
tum distribution, which is a Gaussian centered at zero 
(-P(O) = 0), and whose root mean square is a. This distri- 
bution is assumed to be so narrow {a <C 1), that it only 
populates the first Brillouin zone in a significant way, 

1 / 1 \ 

"^ eXP l"^J 5 "°- (21) 

Note that a ~ 1CP 2 for Bose-Einstein condensates. There- 
fore, each quasimomentum /3 is associated to a plane wave 
centered at k(3, and the CSBW distribution tt^^ is thus 
independent of £, 

Figures [2] and [3] show that the behavior associated to 
the central quasimomentum is visible in the early dy- 
namics. Then, it disappears after a characteristic time 
Td, which decreases as the distribution width a increases. 
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Namely, Fig. [5] shows a numerical simulation of the res- 
onant case k = 4-7T in which this modification manifests 
itself as a transformation of the initial ballistic motion, 
characterized by a quadratic increase of the average ki- 
netic energy (see inset) into a diffusive motion in later 
times t > Td, characterized by a linear increase of the 
average kinetic energy. For the sharper initial distribu- 
tion (in dashed lines), the change in the dynamics oc- 
curs slightly later, and the resulting diffusion is character- 
ized by a higher rate. Fig. [3] shows the anti-resonant case 
k — 2-7T, whose characteristic oscillations of the average 
kinetic energy are damped. The broader is the initial mo- 
mentum distribution, the faster is the damping process. 

Such changes between the early and the asymptotic dy- 
namics has been observed for thermal gases, i.e. incoher- 
ent mixtures of plane waves [281144] . whereas our simula- 
tions are made for Bose- Einstein condensates, i.e. coherent 
superpositions of plane waves. For narrow momentum dis- 
tributions (a <C 1), the dynamics is the same for coherent 
and incoherent superpositions. As the initial momentum 
and quasimomentum distributions are identical, the dif- 
ferent momentum components evolve independently from 
each other, and so are not sensitive to any phase difference 
between them. 

In order to give a simple picture of the dynamics, in 
what follows, we consider a square initial momentum dis- 
tribution 




for \$\ < f 
otherwise. 



(23) 



The parameter A is such that the Gaussian and square 
momentum distributions have the same average kinetic 
energy, that is a — AjyYl. The corresponding CSBW 
distribution n^p is 




4000 



c 3000 
o 

o 

~ 2000 

"l— 

£ 

= 1000 



I I 

I \ ' \ j\ 
' \ ' I 1 \ 
— 1 \ 1 \ 6 \ 

'I | I S / | ' 


i 

— i — i- — 
\; i 


i \ 

i 

i ' ( 


i ' i ' / 1 ; 
■ ; 1 1 i / 1 J 
1 1 i j / \ ij 
+ 1 / i -' 1 9 
i ;,i i J 

1 i y\ \ 1 


I. 5 

i i \ ji 

J J \ i 

1 ' ' 1 
\ i i 

\ ! 

\K i 


i \~ 

i \ 

i \ 
i I 

V. j 1 

1 ' J - 



-0.04 



-0.02 





p 



0.02 



0.04 



Figure 4. The filtering function fp(t) [see Eq. ((25)) ] in full 
lines, and its approximation /^(t) [see Eq. (|27[) [ in dashed lines: 
for (a) t = 8 kicks, and (b) t = 60 kicks. The initial momentum 
distributions, of Gaussian shape (a = 0.0115) and of square 
shape (A = 0.04), are also plotted for comparison, in dotted 
and dashed-dotted lines, respectively. They are multiplied by 
an arbitrary factor to be visible on the graph. 



for \p\ < f 
otherwise. 



(24) 



Applying Eq. (|17p to this particular case, we can write the 
average kinetic energy in the form 



E(t) = E(0) + g 



A/2 



-A/2 



(25) 



which shows a filtering process in momentum space, asso- 
ciated to the function 



sin 2 (ir£ (13 + 1/2) t) 
sin 2 (ni (/3 + 1/2)) 



(26) 



characteristic of diffraction by a grating [23J ■ As fp (t) be- 
comes sharper with time, it better "resolves" the details 
of the initial momentum distribution, especially its finite 
width, which can explain the damping process observed 
on Figs. [Hand El 

An alternative explanation of this damping effect can 
be obtained by considering the behavior of a slice at po- 
sition £ of the initial distribution. After t kicks, the parts 



of the distribution characterized by the extremal values 
of (3 = ±A/2 have been translated, according to (JSJ, of 
ttI (±A + 1) t/2. The filtering process becomes significant 
when an initial slice £ has spanned one potential well (the 
"spatial'' first Brillouin zone), which gives a damping time 
Td = 2/£A ~ (£A)~ ~ (%A)~ , a trend confirmed by our 
numerical simulations. 

Now let us turn to the specific case of Fig. [5] Let us 
approximate the central lobe t < 1 of the filtering func- 
tion f fj (t) by 

f'p(t) =t 2 cos 2 (n(3t), (27) 

which is shown in Figs. |Ua) and (b) for different times 
t = 8 and t = 60. If t < A -1 , the central lobe completely 
encompasses the initial momentum distribution; therefore 
the average kinetic energy is 



E(t) » E(0) 
«E(0) 



K 2 t 2 



AA 



A/2 



d/3 cos 2 (itpt) 

1 -A/2 

K 2 t 2 sin(TTZit)' 



1 



TlAt 



(28) 



M. Lepers, V. Zehnle, J.-C. Garreau: Kicked-rotor quantum resonances in position space. 



7 



which shows, as expected, that the early dynamics is bal- 
listic (E(t) cx t 2 ) 1 with a coefficient K 2 /4, as for a plane 
wave. As t approaches A^ 1 , Eq. (|28p predicts that the dy- 
namics remains ballistic, but that the rate decreases from 
K 2 /A to K 2 /8. This conclusion in not valid for t > A~ x , 
when the edges of np £ fall outside of the central lobe of the 
filtering function. We can then calculate the integral (|25l) 
over the range of the central lobe, that is for |/3| < l/2t: 



E(t) « £7(0) 



KH 1 

~Ja 



l/2t 



l/2t 



dp cos 2 (ir/3t) , (29) 



which yields a diffusive behavior (£7(i) oc t) 



E(t) « E(0) 



8A 



t, 



(30) 



whose rate is inversely proportional to A. This trend is 
also visible on Fig. [2] for Gaussian momentum distribu- 
tions. 

In the anti-resonant case shown in Fig. [3J the energy, 
although damped, oscillates for a time equal to a few A^ 1 , 
and finally "freezes" at K 2 /8. To catch the most impor- 
tant features of the dynamics, we will once again take 
an approximate form for Eq. (|25p . Indeed, we can take 
sin 2 {n£ (/? + 1 /2)} ps 1 over the width of the quasimomen- 
tum distribution A. By expanding sin 2 {tt£ (j3 + 1/2) t}, 
one gets 



E(t) « E(0) + ^ 

- K 2 
« £7(0) + — 



1 

l -A 



l-(-l) 



A/2 
-A/2 

t sin (it At) 



d/3cos{27ri (£ + 1/2)} 
(31 



For t ^> A 1 the oscillations are completely damped and 
the energy tends to its asymptotic value 



E(t -> oo) = £7(0) + 



X 2 



(32) 



Incidentally, if one keeps in Eq. (|3"Tj) a general momen- 
tum distribution irp £, the equation shows that the kinetic 
energy is proportional to the cosine-Fourier series of the 
momentum distribution (see Fig. [3J) . Using the properties 
of Fourier series, we can "invert" Eq. ([31jl to obtain 



MP) 



E 4s « 

s=l 



E(t)) cos ( 2wt 



"A 



(33) 

which yields an approximate expression for the initial dis- 
tribution. Since it appears as a Fourier series, the accuracy 
of this approximation increases when t — > oo. This result 
is valid for narrow momentum distributions that contain 
no ballistic quasimomentum. As ballistic quasimomenta 
are more numerous with increasing k, the phenomenon is 
more likely to be observed for k = 2tt. 



6 Broad initial momentum distribution 

For the sake of completeness, we consider now the trivial 
case of an initial momentum distribution that is a few 
times larger than the Brillouin zone. This case leads to a 
well-known diffusive behavior, 



(34) 



which was demonstrated with thermal gases [23 , 28 , 44 . In 
this section, we use our approach to demonstrate that the 
diffusion can be observed both with coherent and incoher- 
ent superpositions of plane waves. 

The momentum-space wavepacket associated with a 
coherent, say Gaussian, superposition of plane waves is 



-i(n+l3)cj> 



exp 



4ct 2 



(35) 



where a ^> 1. The corresponding CSBW distribution is 
approximately uniform is /3, but strongly localized about 
£ = (f>. Eq. (JSJ shows that, after one kick, the part of the 
CSBW distribution corresponding to quasimomentum /3 
moves to + 2ir£ ((3 + 1/2), which means that the CSBW 
distribution covers a whole potential well. Therefore, the 
average momentum does not change, and the average ki- 
netic energy is increased by J d/3 K 2 sin 2 (£ + 2ir£/3)/2 = 
K 2 /A, hence the diffusive behavior. 

For incoherent mixtures, there is no well-defined phase 
<j) between the different momentum components. By apply- 
ing Eq. (|15p and by supposing that all quasimomenta are 
equally populated, we can write irp t £ s» (2-7r) . The rea- 
soning made for a coherent superposition of plane waves 
can be applied for a given slice £, hence the diffusive be- 
havior characterized by the rate K 2 /4, which is then not 
modified by the average over £. 



7 Conclusion 

In this paper, we have applied the analysis of quantum 
resonances in position space given in [57] to cases that are 
met in the everyday laboratory life. Our interpretation of 
the simple quantum resonances, based on the local mo- 
mentum transferred from the kicks, can also be applied to 
wave functions that are not sharply localized in the "spa- 
tial first Brillouin zone". As shown above, this allows one 
to build intuitive and easily understandable images of the 
phenomena that are often simpler than the corresponding 
explanation in the momentum representation. We have 
applied our approach to both coherent and incoherent su- 
perpositions of plane waves, which stand for Bose- Einstein 
condensates and thermal gases, respectively. We have demon- 
strated that the two types of superpositions induce the 
same dynamical behavior, in the two limiting cases of nar- 
row and broad initial momentum distributions. By con- 
trast, we expect significant differences to appear if the 
width of the distribution is comparable to the first Bril- 
louin zone. 
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The case of high-order quantum resonances has not 
been addressed here. For U — it, the initial wave packet 
reconstructs into two sub-packets separated by half the 
lattice step. When a kick is applied, the momentum shift 
is due to the local potential gradient at the place of the 
sub- packets [27]. However, the weights of the sub-packets 
depend on their phase difference, giving birth to a com- 
plex interfering pattern after a few kicks. The evolution 
of the average kinetic energy reflects this "history", and 
therefore cannot be explained in terms of the simple ar- 
guments developed in this paper. 

An interesting prospect of the present work is its ap- 
plication to the case where interactions between atoms are 
present (e.g. one uses a dense Bose-Einstein condensate). 
A first experimental study of such as system has been 
made by Raizen and co-workers |50| . and it was demon- 
strated in [49] that for a weak interaction strength, the res- 
onant value of k is shifted. The latter work also shows that 
the diffusion coefficient presents sharp edges as the non 
linearity is varied. A clear picture of the physics underly- 
ing these observations is still missing. A kind of nonlinear 
optical Talbot effect has been observed recently in an op- 
tical system [53] . which hints that our approach could be 
generalized to this much more complex case. 
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